Coherent Quantum Dynamics of a Superconducting Flux Qubit 
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We have observed coherent time evolution between two quantum states ol a superconducting flux 
qubit comprising three Josephson junctions in a loop. The superposition of the two states carrying 
opposite macroscopic persistent currents is manipulated by resonant microwave pulses. Readout 
by means of switching-event measurement with an attached superconducting quantum interference 
device revealed quantum-state oscillations with high fidelity. Under strong microwave driving it was 
possible to induce hundreds of coherent oscillations. Pulsed operations on this first sample yielded 
a relaxation time of 900 nanoseconds and a free-induction dephasing time of 20 nanoseconds. These 
results are promising for future solid-state quantum computing. 



It is becoming clear that artificially fabricated solid- 
state devices of macroscopic size may, under certain con- 
ditions, behave as single quantum particles. We report 
on the controlled time-dependent quantum dynamics be- 
tween two states of a micron-size superconducting ring 
containing billions of Cooper pairs (1). From a ground 
state in which all the Cooper pairs circulate in one direc- 
tion, application of resonant microwave pulses can excite 
the system to a state where all pairs move oppositely, 
and make it oscillate coherently between these two states. 
Moreover, multiple pulses can be used to create quantum 
operation sequences. This is of strong fundamental inter- 
est because it allows experimental studies on decoherence 
mechanisms of the quantum behavior of a macroscopic- 
sized object. In addition, it is of great significance in the 
context of quantum computing (2) because these fabri- 
cated structures are attractive for a design that can be 
scaled up to large numbers of quantum bits or qubits (3). 

Superconducting circuits with mesoscopic Josephson 
junctions are expected to behave according to the laws 
of quantum mechanics if they are separated sufficiently 
from external degrees of freedom, thereby reducing the 
decoherence. Quantum oscillations of a superconducting 
two-level system have been observed in the Cooper pair 
box qubit using the charge degree of freedom (4). An 
improved version of the Cooper pair box qubit showed 
that quantum oscillations with a high quality factor could 
be achieved (5). In addition, a qubit based on the phase 
degree of freedom in a Josephson junction was presented, 
consisting of a single, relatively large Josephson junction 
current-biased close to its critical current (6,7). 

Our flux qubit consists of three Josephson junctions 
arranged in a superconducting loop threaded by an exter- 
nally applied magnetic flux near half a superconducting 
flux quantum $ = h/2e [(8); a one-junction fiux-qubit 
is described in (9)\. Varying the flux bias controls the 
energy level separation of this effectively two-level sys- 
tem. At half a flux quantum, the two lowest states are 
symmetric and antisymmetric superpositions of two clas- 



sical states with clockwise and anticlockwise circulating 
currents. As shown by previous microwave spectroscopy 
studies, the qubit can be engineered such that the two 
lowest eigenstates are energetically well separated from 
the higher ones (10). Because the qubit is primarily bi- 
ased by magnetic flux, it is relatively insensitive to the 
charge noise that is abundantly present in circuits of this 
kind. 

The central part of the circuit, fabricated by electron 
beam lithography and shadow evaporation of Al, shows 
the three in-line Josephson junctions together with the 
small loop defining the qubit in which the persistent 
current can flow in two directions, as shown by arrows 
(Fig. 1A). The area of the middle junction of the qubit is 
a = 0.8 times the area of the two outer ones. This ratio, 
together with the charging energy Ec = e 2 /2C and the 
Josephson energy Ej = hlc/^Tre of the outer junctions 
(where Iq and C are their critical current and capaci- 
tance, respectively), determines the qubit energy levels 
(Fig. 2A) as a function of the superconductor phase 7 9 
across the junctions (Fig. IB). Close to j q — ir, the loop 
behaves as a two-level system with an energy separation 
Eiq = Ei — Eq of the eigenstates |0) and |1) described by 
the effective Hamiltonian H = —ea z /2 — Aa x /2, where 
o z%x are the Pauli spin matrices, A is the level repulsion, 
and e w I p $o(j q — tt)/tt (where I p m 2naEj/^ is the 
qubit maximum persistent current) (11). 

The sample is enclosed in a gold-plated copper shield- 
ing box kept at cryogenic temperatures T = 25 mK (k#T 
<C A). The qubit is initialized to the ground state sim- 
ply by allowing it to relax. Coherent control of the qubit 
state is achieved by applying resonant microwave exci- 
tations on the microwave (MW) line (Fig. IB), thereby 
inducing an oscillating magnetic field through the qubit 
loop. The qubit state evolves driven by a time-dependent 
term (— l/2)e mw cos(2TrFt)a z in the Hamiltonian where 
F is the microwave frequency and e mw is the energy- 
modulation amplitude proportional to the microwave am- 
plitude. This dynamic evolution is similar to that of spins 
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Josephson junctions of critical current ~2.2 mA). Evaporat- 
ing Al from two different angles with an oxidation process be- 
tween them gives the small overlapping regions (the Josephson 
junctions). The middle junction of the qubit is 0.8 times the 
area of the other two, and the ratio of qubit/SQUlD areas is 
about 1:3. Arrows indicate the two directions of the persistent 
current in the qubit. The mutual qubit/SQUID inductance 
is M m 9 pH. (B) Schematic of the on-chip circuit; crosses 
represent the Josephson junctions. The SQUID is shunted by 
two capacitors (~5 pF each) to reduce the SQUID plasma fre- 
quency and biased through a resistor (~150 ohms) to avoid 
parasitic resonances in the leads. Symmetry of the circuit 
is introduced to suppress excitation of the SQUID from the 
qubit-control pulses. The MW line provides microwave cur- 
rent bursts inducing oscillating magnetic fields in the qubit 
loop. The current line provides the measuring pulse and 
the voltage line allows the readout of the switching pulse V ou t ■ 
The V ut signal is amplified, and a threshold discriminator 
(dashed line) detects the switching event at room tempera- 
ture. 



in magnetic resonance. When the MW frequency equals 
the energy difference of the qubit, the qubit oscillates 
between the ground state and the excited state. This 
phenomenon is known as Rabi oscillation. The Rabi fre- 
quency depends linearly on the MW amplitude (12-14). 

Readout is performed with an underdamped super- 
conducting quantum interference device (SQUID) with 
a hysteretic current-voltage characteristic in direct con- 
tact with the qubit loop (Fig. 1A). The mutual coupling 
M is relatively large because of the shared kinetic and ge- 
ometric inductances of the joint part enhancing the qubit 
signal. After performing the qubit operation, a bias cur- 
rent pulse lb is applied to the SQUID (15). The I u pulse 
consists of a short current pulse of length ~50 ns followed 
by a trailing plateau of ~500 ns (Fig. IB). During the cur- 
rent pulse, the SQUID either switches to the gap voltage 
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FIG. 2: (A) Calculated energy diagram for the three-junction 
qubit, for Ej/E c = 35, E c = 7.4 GHz and a = 0.8 (11). A-y q 
indicates the phase shift induced by the SQUID bias current. 
(B) Ground-state transition step: The sinusoidal background 
modulation of the SQUID (h g ) is subtracted from the /(, pulse 
amplitude corresponding to 50% switching probability (I sw ) 
and then normalized to I step, the middle value (at the dashed 
line). A sharp peak and dip are induced by a long (1 fjg) MW 
radiation burst at 16 GHz, allowing the symmetry point to be 
found (midpoint of the peak/dip positions, dotted line). Data 
show Isw versus A<& ex t, the deviation in external flux from 
this point. The transition step is displaced from this point 
by A^ q /2n. (C) Frequency of the resonant peaks/dips (dots) 
versus A& ext ; the continuous line is a numerical fit with the 
same parameters as in (A) leading to a value of A = 3.4 GHz, 
whereas the dashed line depicts the case A = 0. 



or stays at zero voltage. The pulse height and length are 
set to optimize the distinction of the switching proba- 
bility between the two qubit states, which couple to the 
SQUID through the associated circulating currents. Be- 
cause the readout electronics has a limited bandwidth of 
^100 kHz, a voltage pulse of 50 ns pulse is too short to be 
detected. For that reason the trailing plateau is added, 
with a current just above the retrapping current of the 
SQUID. The whole shape is adjusted for maximum read- 
out fidelity. The switching probability is obtained by re- 
peating the whole sequence of reequilibration, microwave 
control pulses and readout typically 5000 times. 

When the SQUID bias current is switched on, the cir- 
culating current in the SQUID changes. This circulating 
current, coupled to the qubit through the mutual induc- 
tance, changes the phase bias of the qubit by an esti- 
mated amount 0.017T. Consequently the phase bias at 
which the quantum operations are performed is different 
from the phase bias at readout. This can be very use- 
ful because at the phase bias near n, where the qubit is 
least sensitive to flux noise, the expectation values for 
the qubit circulating current are extremely small. The 
automatic phase bias shift can be used to operate near tt 
and to perform readout at a bias with a good qubit signal 
(11). Care must be taken that the fast shift remains adi- 
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FIG. 3: (A) Rabi oscillations for a resonant frequency F = 
Eio = 6.6 GHz and three different microwave powers ^4 = 0, 
-6 and -12 dBm, where A is the nominal microwave power 
applied at room temperature. The data are well fitted by 
exponentially damped sinusoidal oscillations. The resulting 
decay time is ~150 ns for all powers. (B) Linear dependence 
of the Rabi frequency on the microwave amplitude, expressed 
as 1G" 4 / 20 . The slope is in agreement with estimations based 
on sample design. 



abatic and that the whole sequence is completed within 
the relaxation time. 

The average SQUID switching current I sw versus ap- 
plied flux shows the change of the qubit ground-state cir- 
culating current (Fig. 2B). Here, the pulse amplitude 
is adjusted such that the averaged switching probabil- 
ity is maintained at 50%. A step corresponding to the 
change of qubit circulating current was observed (around 
the dashed line). The relative variation of 2.5% of I sw 
is in agreement with the estimation based on the qubit 
current I p and the qubit-SQUID mutual inductance M . 

The relevant two energy levels of the qubit were first 
examined by spectroscopic means. Before each readout, 
a long microwave pulse (1 fis) at a series of frequencies 
was applied to observe resonant absorption peaks/dips 
each time the qubit energy separation E w — adjusted 
by changing the external flux — coincides with the MW 
frequency F (10). The dots in Fig. 2C are measured 
peak/dip positions, obtained by varying F, whereas the 
continuous line is a numerical fit produced by exact di- 
agonalization (compare Fig. 2A) giving an energy gap 
A w 3.4 GHz. The curves in Fig. 2, B and C, are 
plotted against the change AQ ext in external flux from 
the symmetry position indicated by the dotted line. 
In agreement with our numerical simulations, the step 
(Fig. 2B) is shifted away from the symmetry position 
of the energy spectrum (Fig. 2C) by a phase bias shift 
« 27r(A$ ea;t /$o) « 0.008rr. The step reflects the 
cxtcrnal-flux dependence of the qubit circulating current 
at I sw (after the shift), whereas the spectrum rc- 
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FIG. 4: (A) Ramsey interference: The measured switching 
probability (dots) is plotted against the time between the 
two 7r/2 pulses. The continuous line is a fit by exponentially 
damped oscillations with a decay time of 20 ns. The Ramsey 
interference period of 4.5 ns agrees with the inverse of the 
detuning from resonance, 220 MHz. The resonant frequency 
is 5.71 GHz and microwave power A = dBm. (B) Spin- 
echo experiment: switching probability versus position of the 
7r pulse between two 7r/2 pulses. The period of ~2.3 ns cor- 
responds well to half the inverse of the detuning. The width 
and timing of microwave pulses in the MW line are shown in 
each graph. The readout pulse in the bias line immediately 
follows the last n/2 pulse (see Fig. IB). 



fleets E i at lb = (before the shift) (16). 

Next, we used different MW pulse sequences to induce 
coherent quantum dynamics of the qubit in the time do- 
main. For a given level separation E w , a short resonant 
MW pulse of variable length with frequency F — E w was 
applied. Together with the MW amplitude, the pulse 
length defines the relative occupancy of the ground state 
and the excited state. Corresponding switching proba- 
bility was measured with a fixed-bias current pulse am- 
plitude. We obtained coherent Rabi oscillations of the 
qubit circulating current for a frequency F = 6.6 GHz 
and three different values of the MW power A (Fig. 3). 
The variation in switching probability is around 60%, in- 
dicating that the fidelity in a single readout is of that 
order. By varying A, we verified the linear dependence 
of the Rabi frequency on the MW amplitude, a key signa- 
ture of the Rabi process (Fig. 3B). The oscillation pattern 
can be fitted to a damped sinusoid. For relatively strong 
driving (Rabi period below 10 ns), decay times TR a u up 
to ^150 ns are obtained. This large decay time resulted 
in hundreds of coherent oscillations at large microwave 
power. 

The Rabi scheme also allows the study of the state 
occupancy relaxation. This can be done by applying a 
coherent ir pulse for full rotation of the qubit into the 
excited state and varying the delay time before readout. 
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Experiments performed at F = 5.71 GHz gave an expo- 
nential decay with relaxation time r re i ax « 900 ns. 

As a next step we measured the undriven, free- 
evolution dephasing time by performing a Ramsey in- 
terference experiment (17) as follows. Two n/2 pulses, 
whose length is determined from the Rabi precession pre- 
sented above, are applied to the qubit. The first pulse 
creates a superposition of the |0) and |1) states. If the 
microwave frequency is detuned by SF = E w — F away 
from resonance, the superposition phase increases with 
a rate 2nSF, in the frame rotating with the MW fre- 
quency F. After a varying delay time, we apply a second 
7r/2 pulse to measure the final |0) and |1) state occu- 
pancy via the switching probability. The readout shows 
Ramsey fringes with a period 1/SF, as in Fig. 4A, where 
E w = 5.71 GHz and SF = 220 MHz. The dots repre- 
sent experimental data, whereas the continuous line is an 
exponentially damped sinusoidal fitting curve, yielding a 
free-evolution dephasing time w 20 ns. Note that the 
oscillation period of 4.5 ns agrees well with 1/SF. 

Additional information on the spectral properties of 
the decohering fluctuations can be obtained with a mod- 
ified Ramsey experiment. By inserting a ir pulse between 
the two 7r/2 pulses (Fig. 4B), we obtain a spin-echo pulse 
configuration. The role of the tt pulse is to reverse the 
noise-driven diffusion of the qubit phase at the midpoint 
in time of the free evolution. Dephasing due to fluctu- 
ations of lower frequencies should be cancelled by their 
opposite influence before and after the tt pulse ( 18). Spin- 
echo oscillations (Fig. 4B) are taken under the same con- 
ditions as the Ramsey fringes, but are here recorded as a 
function of the tt pulse position. The period (~2.3 ns) is 
half that of the Ramsey interference. We measured the 
decay of the maximum spin-echo signal (i.e. with the 
tt pulse in the center) versus the delay time between the 
two 7r/2 pulses. The data can be fitted to a half-Gaussian 
(not shown) with a decay time T ec ho ~ 30 ns. 

We conclude that with the present device and setup, 
the dephasing time « 20 ns, as measured with the 
Ramsey pulses, is much shorter than the relaxation time 
Treiax ~ 900 ns. Dephasing is probably caused by a varia- 
tion in time of the qubit energy splitting, due to external 
or internal noise. A likely source is external flux noise, 
which can be reduced in the future. The present qubit 
could not be operated at the symmetry point j q = tt 
where the influence of flux noise is minimal (5), presum- 
ably as the result of an accidentally close SQUID res- 
onance (19). Other possible noise sources are thermal, 
charge, critical current and spin fluctuations. From esti- 
mations of the Johnson noise in the bias circuit (20,21), 
we find a contribution that is several orders of magnitude 
weaker. 

For strong driving, Rabi oscillations persisted for times 
much longer than r^. This constitutes no inconsistency. 
The dependence of the Rabi period on the detuning, due 
to fluctuations of the qubit energy Eiq, is weak when 



the Rabi period is short. The fact that coherence is 
only marginally improved by the tt pulse in the spin- 
echo experiment seems to indicate the presence of noise 
at frequencies beyond 10 MHz. Further analysis and ad- 
ditional measurements are needed. 

These first results on the coherent time evolution of 
a flux qubit are very promising. The already high fi- 
delity of qubit excitation and readout can no doubt be 
improved. Quite likely it is also possible to reduce the de- 
phasing rate. Taken together, these results establish the 
superconducting flux qubit as an attractive candidate for 
solid-state quantum computing. 
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